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Motivated by the experimental search for Bose condensation of quasiparticles in semiconductors, 
the response functions of a weakly interacting Bose gas, with isotropic but non-quadratic dispersion, 
are considered. Non-quadratic dispersion modifies the definition of particle current, and leads to 
modified sum rules for the current-current response function. The effect of these modifications on 
the Berezhinski-Kosterlitz-Thouless transition is discussed. 
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Recently, there has been increasing interest in Bose 
condensation of quasiparticles in solid state systems. 
Examples include indirect excitons in semiconductor 
quantum wells 1 , exciton-polaritons in semiconductor 
microcavities^2*^, quantum hall bilayer excitons^, and 
spin "triplons" in copper compounds 7,8-9 . In many of 
these cases, the composite nature of the quasiparticle 
leads to significant deviations from a quadratic disper- 
sion. Such deviations mean that a current defined by 
J(x) — ip< (x)iVip(x) is no longer correct: such a cur- 
rent is not conserved, and so its correlation functions do 
not obey simple sum rules. Neither can this problem 
be extricated by working in terms of more fundamen- 
tal fields, e.g. the photon/exciton fields for the polariton 
problem, as in such an example the photon current is not 
conserved, the Hamiltonian has terms by which photon 
current is transfered to exciton current and back again. 

There is an obvious solution to this problem: the cor- 
rect definition of current is the Noether current associ- 
ated with invariance of the action under global phase 
rotations. Such a definition automatically leads to a 
conserved current, which for quadratic dispersion is just 
the standard definition. The definition of current and 
its response functions are of particular importance due 
to another common feature of these solid state systems 
in which condensation is sought: they are two dimen- 
sional, and therefore the transition is of the Berezhinski- 
Kosterlitz-Thouless (BKT) classic. Therefore, to find 
the critical conditions at which the transition should oc- 
cur, it is necessary to find the superfluid stiffness, in- 
cluding effects of depletion by density fluctuations. This 
is most naturally achieved by finding the current re- 
sponse functions, and thus separating the current re- 
sponse into normal and superfluid componentsi^ii&ii. 
For the weakly interacting case, one may then pertur- 
batively evaluate the current response functions. Such 
a perturbative evaluation relies on two properties of the 
current response: a sum rule on the longitudinal response 
function (a consequence of using a conserved current)^, 
and an understanding of the effect of vertex corrections 
on the transverse response functionsi£iiLi£. 

The aim of this article is to discuss the correct gen- 
eralisation of response functions for non-quadratic, but 
isotropic, quasi-particle dispersion. Previous work on the 



BKT transition in a model of weakly interacting bosons 
with non-quadratic dispersioni2*2S did not generalise the 
current in this manner. As a result, the current in that 
work was not conserved, and so there is no sum rule relat- 
ing the longitudinal response function to density. As the 
spectrum considered there was quadratic for small mo- 
menta, any formalism which recovers the standard form 
at low densities (i.e. when only low momentum parti- 
cles are excited) will agree. However, at higher den- 
sities, when particles beyond the quadratic dispersion 
contribute to the current response, there are differences 
between the method described here, and the method in 
those previous works, as will be shown below. 

To be precise, consider the following model of a weakly 
interacting Bose gas: 

H = E + f E 4+q^-q^k, 

k k,k',q 

where £k is isotropic, and has a quadratic part as k — > 0, 
but is otherwise general. This Hamiltonian is invariant 
under global phase rotations, and so there is an associ- 
ated Noether current J given by^>: 

Ji{x)= SS iip f (x)- J S r iip(x), (2) 



Sdiip^(x) 



5diip(x) 



where S is the action from the Hamiltonian in eq. 
By definition, this current is conserved, so: 

[H,p(q)] = q- J(q), p(q) = E ^k+q/a^k-q^- ( 3 ) 



In the following, we will be interested in the static re- 
sponse of the system to an applied force that couples to 
such a current, described by the response function: 



Xij(u = 0,q) = 2 / drdJi^^Jji-^O))), (4) 
Jo 

where double angle brackets indicate quantum and ther- 
mal averaging. For an isotropic system, the most general 
form of the response function is: 



Xii(q) = Xr(q) $ 



T J T 



(5) 
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The standard rotating bucket argument still applies in 
dividing the response into a superfluid part that con- 
tributes only to xl and a normal part that contributes 
to both \L and xt- With a quadratic dispersion, the rel- 
evant quantity is p 8 /m = lim g ^o [xl{q) — Xt(q)]- With 
non-quadratic dispersion mass is now q dependent, so the 
identification of p s and m separately is not possible, but 
it is not necessary; the effective vortex action depends 
only on the well defined quantity: 

Xs = lim[xz,(g) - XT(q)}- (6) 

q^Q 

Since the current used is by definition conserved, xl will 
be subject to a sum rule; a generalisation of the sum 
rule that would for quadratic dispersion relate xl to the 
density. Below, this sum rule is evaluated, and thus xl 
and xt are calculated in a perturbative expansion. 

Before evaluating this sum rule, it is first worth stress- 
ing why the above generalisation gives the quantity ap- 
propriate to the BKT transition. The BKT transition 
is associated with the unbinding of vortex pairs. The 
conditions at which the transition occurs are therefore 
described by the effective vortex- vortex interaction, and 
the vortex fugacity. Starting from a microscopic model, 
these quantities both depend on the phase stiffness: the 
coefficient of (V<fi(x)) 2 ~ q 2 <fi 2 in the effective action. It is 
only this quadratic phase dispersion which matters: non- 
quadratic terms in the phase dispersion lead only to short 
range vortex interactions, while the quadratic term leads 
to a logarithmic confining potential. However, the phase 
stiffness is modified by density fluctuations, and the non- 
quadratic dispersion of density fluctuations can modify 
the phase stiffness. Non-quadratic dispersion matters 
because after integrating out density fluctuations, non- 
quadratic dispersion of density fluctuations modifies the 
coefficient of quadratic dispersion of phase fluctuations. 
It is technically easier to evaluate the current response 
functions than to directly integrate out density fluctua- 
tions, and the associated definitions of superfluid stiffness 
are equivalent^. 

The sum rule for xl{q) = QiQjXiji^l) I 'l 2 follows from 
eq. (@J and eq. Q), and the standard procedure, as de- 
scribed for example in refill 

2 f 13 

XL(q) = - / dr((q-J(q,r)q.J(-q,0))) 
<Z Jo 

= ^E e *W[p(q),q-J(q)]k}> (?) 

n 

where one use has been made of the commutation relation 
eq. Writing the commutation relations explicitly in 
terms of the ip \ ip operators, one has: 

q- J(q) = ( £ k+q/2 - Ek-q/2) V' k+q/2 V'k-q/2, 

k 

D>(q)i q • J(q)] = E (£k+ q + e k - q - 2e k ) ^ k+q/2 ^ k _ q/2 - 

k 



In the limit q —> 0, the terms in parentheses are indepen- 
dent of the direction of fc, so we may average over solid 
angles cifi and thus have: 

\hn XL (q) = X>((4V>k)), (8) 

q-> k 

9k = J irH{ ? )■ (9) 

Since dispersion is isotropic, we may write it as e k = 
f(k 2 ), thus g k = (A/d)k 2 f"(k 2 ) + 2/'(fc 2 ), and as ex- 
pected a quadratic dispersion, f(x) — x/2m reduces to 
gk = 1/m. In general, eq. JHJ can be considered as den- 
sity weighted by effective inverse mass of each momentum 
component. The longitudinal response is thus reduced to 
finding an approximation scheme for the occupation of 
each finite k mode, which will be discussed below. 

As yet we have only written explicitly the longitudinal 
component of the current. To find correlations of the 
transverse component, it is convenient to write: 

•Mq) = E *U7*(k + q, k)* k , * k - ( ) , 

(10) 

and from conservation of current, we have: 

gi7i(k + q,k) = cr 3 (e k+q - e k ) • (11) 

Thus, we know the projection of the vector 7^ onto one 
axis; what remains is to find its direction. This follows 
from the definition in eq. J2J), which shows that under the 

interchange "0k V'k ^ ne curren t changes as Ji — > —J,. 
With a little algebra, it can be seen that this directly 
implies 7i(p,q) = (pi + qi)f (p,q)a 3 where /(p,q) is a 
scalar function chosen to satisfy eq. I|ll|). 

From this definition of current, it is now possible to cal- 
culate the current response function perturbatively. As 
in the quadratic dispersion case, the leading order pertur- 
bative calculation relies on properties of the corrections 
to the current vertex 7^ that result from interactions. For 
clarity, the standard argumentiS*^ is summarised here. 
A full calculation of the response would be given by 

lim Xi ,(g)=2 J A^Tr(g(k)r i (k,k)a(k) 7i (k ) k)), 

(12) 

where <?(k) is the Green's function in the Nambu repre- 
sentation indicated in eq. I|10[) . and Ti is the vertex 7^ 
including corrections. At one loop order, vertex correc- 
tions are necessary to make eq. ()12|l satisfy the sum rule, 
eq. |[HJ. However, it can be seen that these required ver- 
tex corrections are of the form shown in Fig. ^ Since 
these involve a vertex where current couples directly to 
the condensate, they involve 7i(q, 0), which, due to the 
previous discussion of the direction of 7^ is proportional 
to qi. Such a correction therefore only changes the lon- 
gitudinal response. Therefore, we may safely evaluate 
the transverse response at one-loop order without such 
corrections. 
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(7»(qiO)VA>. 




r,(k + q,k) = 7i (k + q,k) + 



FIG. 1: Vertex corrections required at one-loop order. 



Having understood why vertex corrections can be ig- 
nored, the perturbative calculation of xt now follows di- 
rectly; writing: 



7i(k, k) = h lim 

q^O 



£ k+q/2 ~~ £ k-q/2 

k~q" 



<j 3 = 2hf'(k 2 ), 



(13) 

with f(k 2 ) — ek as before for an isotropic mass, we thus 
have: 



lim xt (q) = 

q^O d 



d d k 

(27r) d 



k 2 (2f'(k 2 )YTT (g(k)a 3 g(k)a 3 ) . 

(14) 

To complete the evaluation of \t then requires an explicit 
form for the Green's function. The Bogoliubov one-loop 
approximation for the condensed Green's function is: 



S(w,k) 



1 



-iuo + e k + fx —fi 

—fj, iuu + e k + n 



(15) 

where £/. = \J e k (e k + 2fi) is the Bogoliubov quasiparticle 
energy. Thus, eq. (|14|l becomes: 



Tr {Qa 3 Qa 3 ) 



lim xt(q) 

q^O 



- 0J„ 



■> r d d k fc 2 (2f(fc 2 )) 2 

(2n) d 2 



§n' B (&), (16) 



»b(&)(17) 



Finally, to find xl requires evaluation of the average 
occupation of each k mode in eq. JBJ. In evaluating 
eq. JSJ, as the effects of fluctuations in the presence of a 
condensate are required, the condensate depletion due to 
fluctuations must be included in order to derive a con- 



sistent answer"***. This means po 



must 



include fluctuation corrections, determined by consider- 
ing a chemical potential coupled only to k — modes, or 
equivalently by using the Hugenholtz-Pines relation at 
one loop order: 



Po 



S-Ef 2 ((^» + 5(((^» +h " 



Inserting this in eq. JHJ yields the final form: 



(18) 



lim XL (g) = .go 

q^Q 1 2B 



d d k 



(2?r) d 
(2g -g k ) 



9o- 



2e k + p 



9 k- 



£fe + M - £ft , 9oP £fc - £fe - A* 



2& 



2 Cfe( e fc+M) 



(19) 



(To avoid the ultra-violet divergence associated with 
a delta-function interaction, the standard T-matrix 
renormalisationSi has been performed, thus T 2 b is the 
renormalised two-body T-matrix corresponding to the 
bare interaction g). 

In two dimensions, the BKT transition is found, as 
discussed above, by evaluating eq. (fTTjl and eq. (fTT))) at a 
fixed temperature, and finding the value of p which sat- 
isfies liniq^o [xl(q) — Xt(q)] = 7rfc B T/2. The discussion 
up to now has been for a generic isotropic dispersion. For 
illustration, Fig.[2]shows the results of such a calculation 
with a dispersion appropriate to exciton-polaritons: 



2mi 



' k 2 
2m 2 



Oh 



(20) 



with X/mi = l/mx + l/mp, and l/m2 = 1/mx — 1/mp. 
Parameters are chosen close to those of the experiments 
of ref. 4 in CdTe: exciton mass m x — 0.08to c , photon 
mass mp = 2.58 x 10- 5 m e , O fl = 26meV, and T 2B = 
13meV/10 n cm- 2 . 



— Result of this article 

- ■ Method of Kavokin et al. 



O 

o l e +09 




10 100 
Temperature (K) 



FIG. 2: Comparison of calculation of BKT critical density vs 
te mperatu re for the method discussed here and the method in 
ref. ll9l2Ct At low densities, both calculations agree, as non- 
quadratic effects are irrelevant, but where such effects matter, 
their predictions differ. 

In the normal state, the transverse and longitudinal 
response functions should become equal. It is instructive 
to see how the expression for effective mass, weighting 
the density, appears in such a calculation. In the normal 
state, there are no condensate depletion effects to worry 
about, and so: 



Xl(0) 



d d k 

W) 



j9kn B (f(k )), 



(21) 



The one-loop transverse response is as in eq. (11711 . but 
with £k — * £fc = f{k 2 ). To see that they agree, it is conve- 
nient to rewrite eq. (|17|l with a change of integration vari- 
ables. We first introduce x = fc 2 , so d d k = SdX d ^ 2 ~ 1 dx/2, 



4 



with Sd the surface of the d-dimensional hypersphcre, 
and then change integration variable again to f(x), with 
dx = df/f'(x), giving: 



Xt(0) 



_2 r 

dj (2tt 



S d d/ x^ 2 - 1 x(2f'(x)) 2 dn B (f) 



) d 2f'(x) 
S d df ( 1 



d 



2d r 

d2j 



(2ir) d \2f'(x) dx 
d d k 



2 df 
x d / 2 2f'(x) 



(27T)' 



T9kn B {f(k )), 



(22) 



where the second line is integration by parts, and the last 
used d x (2x d / 2 f(x)) = (d/2)x d / 2 - 1 g(x). 

In conclusion, a formalism for calculating the trans- 
verse and longitudinal response functions of a Bose gas 
with arbitrary isotropic dispersion has been presented. 



A sum rule relates the longitudinal response to density 
weighted by effective inverse mass at a given momen- 
tum. Using such a formalism recovers the equivalence 
of transverse and longitudinal responses in the normal 
state. This formalism allows a consistent formulation of 
the critical conditions for the BKT transition in a two- 
dimensional Bose gas. 
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